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ABSTRACT 


k 


he  oxidation  of  benzene  under  fuel-lean  conditions  has  been 

studied  from  1600  to  2300  K  in  a  shock  tube  using  a  stabilized  cw 

CO  laser  to  monitor  CO  production.  The  result  of  the  kinetic 

modeling  of  the  CO  formation  according  to  the  global  mechanism  of 

Fujii  and  Asaba  indicates  that  the  rate  of  CO  formation  in  the 

early  stage  of  oxidation  depends  very  sensitively  on  the  rate- 

limiting  unlmolecular  decomposition  of  benzene* 

. .  c6h6  i  c6h5  +  h 

It* 


The  rate  constant  obtained  for  the  initiation  reaction  using  a 
scheme  consisting  of  a  set  of  25  reactions  is  given  by  the 
following  expression:  ) 

ki  *  iol5.70iO.52  exp( -54^300  ±  2200/Tf" sec“l. 

'‘for  the  pressure  range  of  1.9  -  2.7  atm^  The  result  of  an  RRKM 
calculation  for  the  reaction  indicated  Cut  at  temperatures  above 
2000  K,  ki  becomes  slightly  pressure-dependent.  The  extrapolation 
of  the  individual  data  with  the  RRKM  theory  led  to 

*1*  *  10l5»76±0»50  exp(-58,400  t  2200/T)  sec~l. 


a.;. 


INTRODUCTION 


Th»r«  it  considerable  enrreat  interest  in  the  kinetics  and 
aechanisa  of  bensene  oxidation  because  of  the  increasing  use  of 
aroaaties  as  fuel  components.  Partly  due  to  the  extreaely 
complex  moehanisms  involved,  to  date  there  exist  only  a  few  high 
teaperatare  kinetic  stadias  (appropriate  to  coabastion)  on  the 
oxidation  of  aroaaties.  Early  shoek  tube  stadias  on  beaxeae 
oxidation1"5  V(r«  concerned  largely  with  ignition  delays  which 
eoald  not  be  readily  interpreted  aechaniatieally. 

Only  in  the  nore  recent  aajor  works  by  Pajii  and  Aaaba^**^^ 
asing  a  shock  tabe  and  by  Claasann  and  co-totk»t*H»^  naing  a 
high  toaporataro  adiabatic  tarbalont  flow  roaetor.  n  elonror 
global  pletaro  of  boaxono  oxidation  began  to  oaorgo. 


One  eoaaon  pronineat  featare  foaad'ia  the  works  of  the 
aforcaoatioaed  two  groaps.  as  well  as  in  the  low  toaporataro 


stady  by  Norrish  and  Taylor15  is  that  CO  is  a  aajor  early 
prodaot.  However,  the  aochnnisa  of  CO  forantioa  renaias  largely 


spoealntivo.  To  accoant  for  the  obsorvsd  CO  foraation,  Fajii  and 
Aaaba1®  proposed  the  'global*  aeehaaisa: 

C4H,  ♦  0, _ _  2C0  ♦  C4H4  «■  H  . 

Olassann  and  co-workers11*55  proposed  a  detailed  aeehaaisa  of 
CO  foraation  involving  the  phonoxy  radical.  The  high  phenol 
concentration  obtained  in  their  work  was  proposed  as  evidence  to 
this  stop.  The  detailed  aeehaaisa  seeas  to  be  roaghly  eonsistent 
p  with  Asaba's  global  scheae.15 


Is  the  present  work,  we  need  e  shock  tube  to  studp  the 
oxidation  of  bensene  under  extremelp  fuel-lean  conditions  under 
which  the  fast  initial  CO  foraatiou  would  be  sensitive  to  the 


initiation  step  of  phenjrl  radical  foraatiou.  The  real  tine  CO 
concentration  is  aeasnred  bp  resonance  absorption  using  a 
f r eque acp- s tab il is ed  CV  CO  laser.  Kinetic  modeling  of  onr  CO 
profiles  as  well  as  its  delap  tines,  was  aade  using  a  set  of  23 
reactions,  based  largelp  on  Fuji!  and  Asaba's  global  seheae  to 
obtain  infornation  on  the  initiation  reaction,  C#H4 — >  C,I,  +  H, 
an  iaportant  and  poorlp  determined  reaction  in  the  bensene 
oxidation  spstea. 

EXPERIMENTAL 

The  details  of  the  shock  tube,  the  f requeaep— stabil lsed  CT 
CO  probe  laser,  CO  absorption  calibration,  and  other  similar 
experiaental  procedures  were  given  in  frPevious  publications.1**15 
la  this  work,  onlp  reflected  shocks  were  used.  All  calculations 
for  shocked  gas  properties  such  as  temperature  and  pressure  were 
made  bp  using  the  NASA/LBWIS  equilibrium  program. **  For  all 
three  reaction  mixtures  with  C(Ha:Oa:Ar  ratios  of  1:200:1799, 
1:214:3011,  and  1:218:6447,  the  2  1  p(10)  transition  of  the  CO 


laser  was  used  to  measure  the  CO  product  concentration  profiles 
bp  resonance  absorption.  CO  threshold  times  were  measured  in 


other  experiments  using  mixtures  with  ratios  of  1:2S0:2070  and 
1:263:2439.  The  laser  output  was  continuouslp  stabilised  to  line 
center  bp  the  use  of  a  loek-ia  s t ab il ixer . ** • *  * Sines 

substantial  concentrations  of  0,  were  used  in  the  CaHa-0a-Ar 

* 

reaction  mixtures,  the  effective  room  temperature  collision 


half-width  r*  tor  CO.  whi< 


broads  slot  ia 


tho  CO  absorption  liaawidth.  was  takaa  to 

T*  “  YC0-0,  X0,  *  YCO-ArXAr 

whara  Xq#  and  X^r  ara  tha  aol«f notions  of  0,  and  Ar. 

Tco-o,  “  °*0452  c“_1  •*■“*  <*•*•  »■<*  r^o-Kt  “  °‘0522 

ca"1  its"1  (Ref.  14)  ara  tha  room  temperature  collision 
half-widths  of  CO  with  0S  and  Ar  as  collision  partners. 


respectively.  The 


lion  of  the  absorptii 


absolute  CO  concentrations  was  achieved  by  the  use  of  the 
calibrated  gain  equation  in  Ref.  14. 


Ve  have  also  examined  tha  possibility  of  UV  omission  by 
electronioally  excited  spooies  snoh  as  C4H(*  and  C4H(0H*  in  the 


2S0  nm  region.  Ernissi* 


ide  by  using  a  1P28 


photomultiplier  with  a  2S4  am  band  filter  which  has  a  FVHM  band 


width  of  10  nm.  Thi 


ilssii 


indeed 


will  bo  shown 


and  discussed  below. 


Certified  A.C.S.  benseae  from  Fisher  Seientifie  was  purified 
by  first  partially  crystalising  it  at  about  279C  in  a  previously 


finger.  Then  the  liquid-solid  mixture  was  pumped 


for  several  minutes  to  preferentially  remove  any  of  tho  more 


volatile  impurities.  The  ail 


tletely  thawed 


and  tho  same  procedure  was  repeated  two  more  times.  0C  analysis 

of  tho  bonsono  samples  before  and  after  purification  showed  a 

large  reduetion  in  the  impurities  and  the  purified  sample  showed 

* 

an  impurity  level  of  only  <  0.002%. 


KINETIC  MODELING  AND  RESULTS 


Reflected  shock  experiments  have  been  carried  out  in  the 
temperature  range  from  1600  to  2313  K  and  pressure  range  from  1.9 
to  2.7  atm.  for  five  CgHg/02/Ar  mixtures  with  CgHg  molefraction: 
0.015%.  0.031%,  0.037%,  0.043%,  and  0.050%  and  the  02/CgHg  ratio 
of  200,  214,  250,  263  and  218. 

Typical  CO  absorption  and  254  nm  emission  traces  are  shown 
in  Fig.  1(a)  and  1(b)  respectively.  It  should  be  mentioned  here 
that  the  delay  time  in  the  emission  correlates  closely  to  that  in 
the  CO  absorption  trace.  Typical  CO  concentration  time  profies  at 
2275,  2076,  and  1743  K  are  shown  by  the  circles  in  Figs.  2(a), 

2(b)  and  2(c)  respectively.  A  conspicuous  feature  of  these 
profiles  is  the  relatively  long  induction  (or  appearance)  times, 
followed  by  rapid  production  of  CO  with  maxima  approaching  the 
total  carbon  limit.  This  suggests  a  high  activation  energy 
initiation  step,  followed  by  fast  radical  reactions  which  rapidly 
convert  the  hydrocarbon  species  (l.e.,  CgHg,  CgHg,  C2  H2 ,  CH2* 
etc.)  to  CO  and  subsequently  to  CO2. 

To  Interpret  our  CO  production  results,  we  carried  out 
kinetic  modeling  using  a  set  of  25  reactions  as  shown  in  Table  X, 
based  largely  on  the  global  mechanism  of  Fujii  and  Asaba, 
established  according  to  their  data  obtained  from  300  nm  UV 
absorption,  5  pm  CO  emission  as  well  as  end  product  analysis  (in 
single-pulsed  shock  experiments).10  We  have  neglected  the  CgHg  + 
°2  *♦  CgHg  +  KO2  reaction  in  our  present  scheme  because 
according  to  the  rate  constant  estimated  from  tratisition-state- 

\ 
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theory  calculations  by  Lin  and  Tevault31,  its  total  reaction  rate, 
even  at  out  O2  concentrations,  at  1600  K  (the  lower  temperature 
limit  in  our  experiment),  is  104  times  less  than  that  of  benzene 
decomposition  (Reaction  >1)).  The  contribution  of  this 
blomolecular  initiation  reaction  was  found  to  be  less  than  5%  of 
the  total  initiation  rate  in  the  worst  case  employed  in  this  work 
when  Fujii  and  Asaba's  unusually  large  rate  constant 

(k  *  6.3  x  1013  e"30'000/T,  cc/mole-sec)1-0 
was  used.  Their  value,  however,  was  established  on  the  basis  of 
the  low  c-H  bond  dissociation  energy  in  benzene,  427  kJ/mole 
(which  is  to  be  compared  with  the  presently  accepted  value,  464  ♦ 

8  kJ/mole  as  will  be  discussed  later). 

The  reactions  C5H5  ->  C2H2  +  C4H3  and  CeHg  +  C5H5 
biphenyl  +  H  were  also  excluded  in  our  modeling  in  view  of  our 
extremely  O2  -  rich  conditions.  The  inclusion  of  these  reactions 
was  found  to  have  little  effect  on  CO  product  profiles  according 
to  the  result  of  out  test  modeling.  The  following  acetylene- 
related  reactions  were  added  to  the  scheme:  C2H2  +  o  *♦  HC2O  +  H, 
C2H2  +  0  CH2  +  CO,  CH2  ♦  O2  ■»  HCO  ♦  OH,  and  HC2O  +  O2  ■»  2CO  + 

r 

OH.  We  have  adopted  Fujii  and  Asaba's  fast  "global"  reaction  for 
CO  production  from  the  C5H5  +  O2  -»  2CO  +  C4H4  +  H  reaction,  as 
this  reaction  was  essential  for  the  successful  modeling  of  our  CO 
profiles,  as  revealed  by  our  sensitivity  analysis  to  be  discussed 
later.  Zn  the  modified  mechanism,  the  only  initiation  step  is 
Reaction  (1),  the  unlmolecular  decomposition  of  benzene.  Zn 
modeling  our  CO  profiles,  it  was  found  that  good  fits  could  be 
obtained  by  varying  only  k^.  Typical  modeling  results  are  shown 


5 


by  the  solid  curves  In  Fig.  2. 


The  quality  of  our  modeling  of  the  CO  formation  profiles  is 
also  indicated  by  the  close  agreement  of  the  experimental  and 
modeled  threshold  (xQ)  as  well  as  the  time  required  to  reach 
half  of  the  maximum  CO  yield  (tj)#  as  shown  in  Fig.  3. 


The  extracted  k^  values  for  benzene  decomposition  obtained 
from  modeling  the  CO  profiles  are  presented  in  the  Arrhenius  plot 
in  Fig.  4.  Also  plotted  in  the  figure  are  k^  values  obtained  from 
modeling  CO  threshold  time  xQ  from  eight  experiments  in  which  only 
threshold  time^nformation  was  available.  A  least  squares  fit  of 

j 

these  points  yields  the  following  expression: 

ki  «  iq15. 70±0.52  exp(-54300  t  2200/T)  sec**. 


Also  included  in  Fig. 4  are  Fujil  and  Asaba's  data,10  covering 
temperatures  from  1200  to  1900  K.  The  agreement  between  the  two 
sets  of  data  is  good  considering  the  drastically  different 
conditions  employed  in  these  experiments.  In  order  to  test  the 
validity  of  our  ki  obtained  from  the  modeling,  a  crude 
sensitivity  analysis  has  been  carried  out  at  2076  K  (run 
corresponding  to  Fig.  2(b))  for  IS  reactions  in  Table  I  which 
include  the  key  reactions  involving  carbon-containing  species. 
Each  of  the  rate  constants  was  varied  by  a  factor  of  two  up  or 
down  and  the  percent  change  in  CO  concentration  at  time  was 
noted.  The  sensitivity  test  results  are  given  in  the  right-hand 
column  in  Table  I.  ^ 
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We  have  carried  out  an  RRKM  calculation  to  check  the  effect 
of  pressure  on  k}  over  the  whole  range  of  temperature  (1600-  2300 


K)  and  pressure  (1.9  -  2.7  atm)  studied  based  on  Troe's  weak 
collision  assumption. 32  It  was  found  that  at  temperatures  above 
2000  K,  ki  becomes  slightly  pressure-dependent.  The  extrapolation 
of  all  data  points  using  an  average  energy  transfer  step-size  of 
l.S  kcal/mole  for  the  oxygen-rich  mixtures  employed  led  to  the 
value  of  kj*  «  10l6*76  *  0.50  exp(-58,400  ♦  2200/T)  sec-1. 

DISCUSSION 

Although  the  mechanism  of  Fujii  and  Asaba8"10  established 
under  fuel  rich  and  lower  temperature  conditions  was  globally 
crude,  it  was  found  to  be  sufficient  to  account  for  our  CO 
production  profiles  measured  under  much  different,  fuel  lean  and 
higher  temperature,  conditions.  The  fitted  rate  constant  for  the 
rate-limiting  initiation  reaction  •  • 

c6h6  -♦  c6H5  ♦  H 

is  reasonable  in  terms  of  the  values  of  both  A-f actor  (6  x  1016 
sec"2)  and  activation  energy  (485  t  18  kJ/mole).  The  latter  gives 
rise  to  the  energy  barrier  at  OK,  Eq  *  463  kJ/mole,  which  is  in 
close  agreement  with  the  recently  accepted  bond  dissociation 
energy,  D(CgH5  -  H)  *  464  t  8  kJ/mole.33 

This  global  mechanism  can  account  not  only  for  the  observed 
CO  production  threshold  and  half -maximum  times,  but  also 
qualitatively  for  the  observed  profiles  of  the  UV  emission  at 
254  nm  as  shown  in  Fig.  5.  This  emission  derives  most  likely  from 

r 

the  electronically  excited  phenol  at  the  state,  which  lies 

\ 
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approximately  443  kJ/mole  above  the  ground  (Sq)  state.34'35  Since 
no  such  emission  was  detected  in  the  shock  heating  of  CgHsOH  or 
C5H5  in  Ar  under  similar  conditions,  the  following  chemical 
pumping  processes  are  energetically  most  likely  candidates  for  the 
production  of  the  excited  CgHsOH: 

0  +  CgHg  -»  CgHsOH*t  #  AH0  =  -429  kJ/mole 
C6H5  +  OH  ■+  C6H50H*t,  AH0  =  -464  kJ/mole 
where  and  "t"  represent  electronic  and  vibrational  excitation 
respectively.  The  energies  available  for  the  electronic 
excitation  at  2000  K  are  approximately  Eav  =  “AH0  +  3RT  = 

480  kJ/mole  in  the  O  +  CgHg  reaction  and  Eav  *  “AH0  +  4  RT 
=  S30  kJ/mole  in  the  CgHs  +  OH  reaction.  Although  both 
reactions  are  not  Included  in  the  global  scheme,  the  profiles 
calculated  from  X  «  (OJ(CgHgl  and  f  C5H5 J [OH]  are  in  qualitative 
agreement  with  the  observed  UV  emission  profiles  as  shown  in 
Fig.  5.  This  finding  is  another  indication  that  the  global 
scheme  of  Fuji!  and  Asaba  is  adequate  for  the  description  of  the 
overall  oxidation  kinetics. 

The  apparent  usefulness  of  the  global  scheme  in  describing 
our  present  findings  lies  perhaps  in  the  correctness  of  the  value 
of  the  rate  constant  assigned  by  Fujil  and  Asaba  for  the  overall 
chain  propagation  step: 

c6h5  +  °2  4  2C0  +  c4h4  +  H» 

This  abridged  step,  as  has  been  discussed  by  Glassman  and 
coworkers,3-2  could  represent  the  following  plausible  sequence 

* 

\ 
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of  events:  C5H5  ♦  <>2  ■*  0  ♦  C5H5O  (0  +  C5H5)  ♦  CO  -*  CO  ♦  C5H5O 
-*  2C0  +  C4H4  ♦  H.  On  the  other  hand,  it  is  also  likely  that 


the  elementary  reaction. 


CgHj  +  02**  CgHgOO* 


01 


■*  2C0  +  C4H4  +  H 

may  occur  at  combustion  temperatures. 


At  present  we  cannot  differentiate  one  electronic  excitation 
process  from  the  other.  He  believe,  however,  that  the  o  +  CgHg 
reaction  is  the  key  process  responsible  for  the  uv  chemilumine¬ 
scence  because  a  direct  test  for  the  reaction  using  N2O  as  the 
0  atom  source  revealed  an  immediate  appearance  of  the  254  nm 
emission  with  no  delay  at  1700  K,  at  which  the  dissociation  of 
CgHg  is  slow.  Additionally,  the  concentration  of  CgHs  may 
be  in  reality  not  as  high  as  that  predicted  by  the  global  scheme 
employed  here.  It  is  quite  possible  that  the  key  chain-carrying 
reactions  involving  the  reaction  CgHg  with  O36  and  OH37  occur 
mainly  by  the  replacement  processes  via  addition,  rather  than  by 
abstraction  yielding  CgHs,  even  at  these  high  temperatures:  viz., 

0  ♦  CgHg  CgHsOH*t  C6H50  ♦  H 
AH°  *  -66.5  kJ/mole 
OH  CgHg  CgHgOH*  CgH5OH  +  H 
AH® -0.8  kJ/mole. 

The  occurrence  of  these  reactions  can  account  for  the  production 
of  large  amounts  of  CgHsOH  and  its  decomposition  products,  CO, 

CgHg ,  C4H6#  etc.,  as  recently  reported  by  Glassman  and 
coworkers. 12 
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Obviously,  the  global  scheme  of  Fuji!  and  Asaba  Is  mechanis¬ 


Y 


tically  too  crude  to  account  for  the  production  of  Cg  and  C5 
Intermediates.  Any  realistic  mechanism  proposed  for  CgHg 
oxdidation  therefore  should  Include  these  key  Intermediates  and 
account  for  their  formation  and  disappearance.  To  achieve  this 
goal,  evidently,  much  more  kinetic  data  are  needed,  particularly 
with  regard  to  various  key  elementary  reactions  involving 
C6(C6Hg,  CgHgOH,  CgH5),  C5(C5H5 ,  C5Hg),  C4(C4H4,  C4Hg)  and  C3HX 
with  0,  H  and  OH  over  a  broad  range  of  temperature  and  pressure. 
Almost  all  of  these  reactions  are  expected  to  occur  via  stable 
adducts  which  may  readily  branch  into  diverse  products  at  high 
temperatures. ~  Further  study  is  underway  at  this  laboratory  for 
the  elucidation  of  the  mechanisms  of  these  processes. 

CONCLUDING  REMARKS 

In  this  work  we  have  studied  the  formation  of  CO  in  the 
oxidation  of  benzene  in  shock  waves  under  fuel  lean  conditions 
employing  the  CO  laser  resonance  absorption  technique.  Kinetic 
modeling  of  CO  production  profiles  based  largely  on  Fujll  and 
Asaba' s  global  mechanism  gave  rise  to  the  extrapolated 
high-pressure  rate  constant  for  the  key  initiation  reaction, 

CgHg  CgH5  +  H: 

ki*  *  6  x  10*6  e-58,400/T  sec”*. 

These  Arrhenius  parameters  appear  to  be  reasonable  in  comparison 
with  the  extrapolated  values  for  toluene  decomposition  under  high 
pressure  conditions. 38 


\ 
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Although  the  reliability  of  k^  thus  derived  depends  on 
the  validity  of  the  mechanism  assumed,  the  "reasonableness"  of  its 
parameters  and  its  ability  to  simulate  CO  production  profiles  over 
the  broad  range  of  conditions  employed  (as  shown  in  Figures  2  and 
3)  seem  reassuring.  The  ultimate  checking  of  this  rate  constant 
would  have  to  rely  on  the  direct  diagnostics  of  either  H  or  C5H5 
formed  in  the  reaction. 
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Reaction*  and  Rate  Conatants  Used  in  the  Modeling 
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Figure  Captious 


Fig.  1.  -  Typical  CO  absorption  trace  (a)  and  2S4  aa  emission 
trace  (b)  at  1792  C  and  2.37  atn  for  a  1:230:2070 
(CcH(:0s:Ar)  mixture.  The  bottom  trace  in  each  panel  is  tba 
pressure  trace. 

Fig.  2.  -  Typical  CO  eonee at ra t ion- t ime  profiles  for  a  1:200:1799 
(C4H<:Oa:Ar)  mixture  at  (a)  2273  K.  2.00  atm.,  (b)  2073  I, 

2.10  atm.,  and  (e)  1743  C,  2.4  atm.  Circles-are - ' 

experimental  data  and  solid  eurves  are  modeled  fits  using 
tbe  scheme  in  Table  I.  The  total  earbon  eontent  in  eaeh 
experiment  is  indicated  by  the  line  denoted  by  Ct#t. 

Fig.  3 .  -  Normalised  CO  threshold  (t#,  eireles)  and  half  maximum 
(t,y,.  triangles)  times  for  a  1:200:1799  (C(Bc:0a:Ar) 
mixture.  Open  symbols  are  experimental  points’  and  filled 
symbols  are  results  from  modeling.  See  text. 

Fig.  4.  -  Arrhenius  plot  of  ka  extraeted  from  modeling  of 
CO  profiles  for  C(Hc:0t:Ar  mixtures  (0-  1:200:1799. 

Q  -  1:214:3011*  A  -  1:218:6447)  and  by  modeling  of 
CO  threshold  times  for  mixtures  (0  -  1:230:2070. 

A  -  1:263:2439).  Result  of  a  least  squares  analysis  is 
shown  by  the  solid  line.  The  dashed  curve  is  Fuji!  and 
Asaba * s  results  walld  from  1200  to  1900  K.  k*  given  by 
the  dotted  line  is  the  extrapolated  high-pressure  rate 

constant  (see  text).  * 
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